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We consider an index-guiding silica photonic crystal fiber with a triangular air-hole structure and a poled 
quadratic nonlinearity. By tuning the pitch and the relative hole size, second-harmonic generation with zero 
group-velocity mismatch is found for any fundamental wavelength above 780 nm. The phase-velocity mismatch 
has a lower limit with coherence lengths in the micron range. The dimensionless nonlinear parameter is inversely 
proportional to the pitch and proportional to the relative hole size. Selected cases show bandwidths suitable 
for 20 fs pulse-conversion with conversion efficiencies as high as 25%/mW. © 2008 Optical Society of America 
OCIS codes: 060.2280, 060.2400, 060.4370, 320.7110 



Relying on quadratic nonlinear it ies, second-harmonic 
generation (SHG) is widely used for efficient wavelength 
conversion devices in order to extend the spectral range 
of laser sources and to do all-optical wavelength multi- 
plexing. Efficient conversion from the the fundamental to 
the second-harmonic (SH) mode requires a small phase 
mismatch between them. Phase matching to the lowest 
order is typically achieved through a quasi-phase match- 
ing (QPM) technique,^ whereby the group- velocity mis- 
match (GVM) sets the limits to device length and band- 
width for pulsed SHG. In conventional fibers, SHG with 
near-zero GVM was found for restricted wavelengths,^ 
while zero GVM was predicted using mode-matching.^ 
For bulk media zero GVM was found for restricted wave- 
lengths by spectrally noncritical phase matching,^ and 
by combining non-collinear QPM with a pulse- front tilt.^ 
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Fig. 1. Sketch of a triangular structured index-guiding 
PCF with pitch A and air-hole diameter d. 

Here we investigate efficient pulsed SHG in a poled 
silica photonic crystal fiber (PCF), having a standard 
index-guiding triangular design with a single rod defect 
in the center (see Fig. The main design parameters 
of the PCF are the pitch A and the relative hole size 
D = d/A. The nonlinearity is induced, e.g., by thermal 
poling as has recently been demonstrated in PCFs.^ 

When the fundamental is assumed undepleted 
and continuous wave, the SHG efficiency is^ r] ex 
PitZg^/pSinc^(A/3/F/2)/Aovi. The contributions to rj of 
nonlinear nature are the fundamental power Pi , the non- 
linear coefficient deff and the effective mode overlap area 
Aov\' Efficient SHG therefore requires a strong mode con- 
finement within a small core and a small effective mode 
overlap area, as well as a strong nonlinear material re- 
sponse. The linear contributions to rj are the fiber length 



/f and the phase mismatch Ap = 2(3i — (^2 between the 
modes {Pj are the mode propagation constants). The 
coherence length of the phase mismatch to the lowest 
order gives the range over which power is exchanged ef- 
ficiently to the SH, but a QPM method can compensate 
for this. The SHG bandwidth is determined by the width 
of sinc^(A/3/F/2), and since ultra-short pulses have broad 
spectra, efficient conversion requires a large bandwidth 
(small Api-p)- Reducing 1^ does the trick, but this gives 
a poor efficiency. The ffist order contribution to Af] is 
GVM, giving a temporal walk-off length Iw between the 
modes, so we must take /f < Iw- When GVM is zero, 
Iw ^ oo, and 2. order dispersion takes over, giving a 
larger bandwidth in a longer device. 

In this Letter, we tune the phase- matching properties 
of SHG by exploiting the flexibility that PCFs offer in 
designing the dispersion properties.^ Previous investiga- 
tions^ of SHG in PCFs considered the scalar case and 
found large bandwidths and strong modal overlaps for 
selected parameter values. Instead, we perform a detailed 
vectorial analysis over a continuous parameter space, 
and show zero GVM for any fundamental wavelength 
Ai > 780 nm by merely adjusting the pitch and the rela- 
tive hole size. This is a much simpler way of removing 
GVM compared to previous methods, and promises 
very large bandwidths due to its flexibility. The nonlin- 
ear properties of the SHG fiber design are also discussed. 

We first focus on the linear dispersion. A fiber mode 
can be described by an effective index n = c/'Uph, ie., 
the ratio of the speed of light c to the phase veloc- 
ity of the mode Vp^ = u;//3. The dispersive charac- 
ter of p gives a phase-velocity mismatch between the 
modes with frequencies uoi (fundamental) and UJ2 = 2cc;i 
(SH), which we classify through the index mismatch 
An = c[l/vph(cJi) - l/vph(^2)] = c[pi/iJi - ^2/^2]' For 
SHG the index mismatch is related to the lowest order 
phase mismatch as A/3 = 47rAn/Ai. The mode group ve- 
locity is instead defined as l/'^g = 9^; = giving a 
GVM (walk-off) parameter di2 = [l/vgiuoi) — l/vg{uo2)]' 

We calculated the dispersion with the MIT Photonic- 
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Bands (MPB) package.^ Each unit cell contained Uq = 
32^ grid points, and the super cell contained ng^ = 5^ 
unit cells. The fundamental mode frequency and group 
velocity were first calculated, followed by iterative calcu- 
lations of the SH until \u;2 — 2c(;i| < 10~^. Subsequently, 
a perturbative approach^^ was used to introduce chro- 
matic dispersion, allowing us to calculate data once over 
a large {D^Pi) parameter space for A unity, and pertur- 
batively calculate the changes as A was varied. 
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Fig. 2. Color maps of (a) di2 and (b) |An| as function of 
D and Ai, keeping the pitch A = 1.6 fim fixed. The solid 
contour in (a) shows di2 = and it is repeated in (b). 

Figure [21 shows the GVM and index mismatch in 
the (Z),Ai) parameter space, keeping the pitch fixed at 
A = 1.6 jam. Along the the solid contour di2 = 0: thus, 
zero GVM is possible for any Ai > 1 fim by choosing 
a proper D. Fig. [S^a) underlines that this is a general 
trend: there the zero-GVM contour is shown for selected 
pitches, and we found di2 = can be achieved for any 
Ai > 0.78 jam. Note that in Fig. Efb) |An| > 0, even 
with non-zero GVM, so the index mismatch cannot be 
zero. This holds also for other A values, so efficient SHG 
will require additional phase-matching such as QPM. 

Figure inijb) shows the values of the index mismatch as 
the zero-GVM contour is traversed. For A = 0.70, 1.0 
and 1.6 jam a cusp appears around Ai :^ A, after which 
I An I increases with Ai. This is because the fundamental 
mode is no longer well-confined in the core while the 
SH, having a smaller wavelength, is still well confined. 
Conversely, for the considered wavelengths the modes 
are always well confined for A = 3.5, 5.0 /im, explaining 
why a small | An| is observed there. The coherence length 
^coh = 7r/A/3 = Ai/(4An) is shown for completeness in 
Fig-IStc), giving typical values in the micron range. 

Let us focus on the telecom, Nd:YAG and Ti: Sapphire 
operating wavelengths (Ai = 1.55 /im, 1.06 /im and 
0.80 /im, respectively.) In Fig. ^a) we then show the 
D-Yslue required to get zero GVM as A is changed. For 
Ai = 0.80 /im zero GVM requires very large I^- values, 
e.g., D = 0.96 for A = 0.70. For such I)- values devi- 
ations from the ideal circular holes must be expected, 
which might infiuence the results presented here. How- 
ever, to our knowledge it is the first time that zero GVM 
has been demonstrated for SHG in any material for this 
wavelength. For Ai = 1.06 /im the lowest required D- 
values are in a range where the ideal round holes should 
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Fig. 3. Extracted zero GVM contours showing (a) Ai, (b) 
I An I and (c) /coh- The curves in (b)-(c) for A = 1.6 /im 
abruptly stop because Ai G [0.5, 2] /im. The arrows indi- 
cate the contour-direction for increasing Ai. 



be preserved. The curves stop for larger A because it is 
no longer possible to get di2 = (it would require D > 1, 
which is unphysical.) For Ai = 1.55 /im a wide range of 
possibilities are offered, but it is advantageous to have 
A < 2 /im because D increases, which leads to higher in- 
tensities due to stronger mode confinement in a smaller 
core diameter dc = A(2 — D). In Fig. ^ we calculate the 
SHG bandwidth AA of the sinc^(A/?/F/2) term by ex- 
pandingi'^ A/3(Ai + AA) = 5];^(m!)-i AA^S^A;^, 
assuming that a QPM grating compensates the m = 
term (as is routinely done in conventional fibers^). Since 
di2 = 0, the 2. order dispersion dominates yielding very 
large bandwidths. Moreover, because di2 = the band- 
width scales as AA ex /p^^^ (instead of AA ex /p^ for 
di2 7^ 0), and thus a longer device can be created with- 
out loosing too much bandwidth. Note also in Fig. ISJa) 
the turn of the zero-GVM contour around D = 0.43 and 
Ai = 1.6 /im, implying that the m = 2 term vanishes, 
giving an increasing bandwidth as observed in Fig.^J^b). 
It is still unclear to what extend the exact location of this 
turning point is infiuenced by parameter uncertainties. 

Using the reductive perturbation method, the di- 
mensionless nonlinear equations for SHG are 
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{dz — iDid^ )ui = iaulu2e 
{dz - di2dt - iD2d^)u2 = ^a^i?/2e'^^"^^ 



(1) 
(2) 



We have assumed a weak nonlinearity [so Ej(r) = 
Aj{z, t)ej{x.)e'^^^^^~^^^'^ +c.c., where x = (x, ?/),] uniform 
over the silica part of the super cell, and weak transverse 
variations in the refractive index so Vx (VxE) = — V^E. 
The "retarded" coordinate z [traveling with velocity 
Vg{uJi)] is normalized to /p, t to the input pulse length 
r, di2 = di2lF/r and Dj = Ip/ {2T'^)d^Pj . We define 
Aj = Uj{2A^lF/NjnjajCrY^^ so Nj{z) = J dt\uj{z,t)\^ 



2 



^ < 



1.0- 
0.9 
0.8 
0,7 
0.6 
0.5 
0,4 

0.3 
180 
160 

Mo- 
no- 

100 
80 
60 
40 



: (a) 



Fundamental wavelength 

X=0.80pm 

^-1.06 urn 

J. =1.55 urn 



pitch A ((jm) 



Fig. 4. Zero GVM contours vs. A with fixed Ais, showing 
(a) D as well as (b) AA for a /p = 10 cm device. 
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lar effect gives the cusp in |An| in Fig. [H^b).] For large 
a decent fundamental mode confinement is observed 
even for A < Ai [compare Fig. EfS) with (1),] giving a 
shift in the peak towards smaller A. Finally, cr gets larger 
when Ai becomes smaller because if d and A are fixed, 
the light is better confined for smaller A. Tabled shows 
choices for efficient SHG with zero GVM using r = 1 ps 
and /p = 10 cm. The large bandwidths imply that pulses 
as short as rum = 21 fs can be converted. The relative 
SHG efficiencies rj/Pi are as high as 25 %/mW. 

To conclude, by tuning the pitch and relative hole size 
in a standard index-guiding silica PCF, we showed that 
SHG with zero GVM is possible for any Ai > 780 nm. 
This is a new and simple way to remove GVM, which 
in addition has conversion bandwidths suitable for down 
to 20 fs pulse-conversion. The SHG nonlinear parameter 
was inversely proportional to the pitch and proportional 
to the relative hole size, due to smaller mode overlap 
areas for lower pitches and larger relative hole sizes, and 
up to 25%/mW conversion efficiencies was found. 

Support from The Danish Natural Science Research 
Council (FNU, grant no. 21-04-0506) is acknowledged. 
M. Bache's e-mail address is bache@com.dtu.dk. 

References 



Fig. 5. Double- log plots of cr vs. A for = 10 cm and 
r = 1 ps. (l)-(4) show energy distributions of the FH 
modes from MPB. nsc = 9 to ensure localized modes. 



is the photon number in the jth wave. The dimensionless 
nonlinear parameter is <j = pIy{'^Tiuj'Iuj2/'i^\'i^2^q(?tY^'^ 
where p = (afas)-^/^! / ^ixe^ (x) • x^'^x) : e2(x)et(x)|, 
X^^^(x) is the Fourier transform of the quadratic nonlin- 
ear dielectric tensor, ej(x) 



3j{x.){soA^lF/hu;jNj) 



1/2 

the dimensionless transverse modes from MPB and aj = 
/ (ix|ej (x) p the mode areas. Integrating over z in Eq. Q 
gives T] = Pip'^llsmc^{A(3lF/2)2ujl/nln2eoc^. Thus, 
is the vectorial equivalent of dlf^/Aovi of Ref. 2. 

The x^^^ tensor is assumed to have the non-zero ele- 
ments xfj\ = Xji] = Xuj = Xui /^^ where i is the main 
direction of the poling voltage, and j is either of the 
2 remaining directions. The present problem has two 
degenerate solutions that are x- and ^-polarized, respec- 
tively, so it will suffice to consider i = x. We used a 

(2) 

realistic value XxJx = 1 pm/V. The nonlinear parameter 
a in Fig. is calculated by fixing Ai and D (in contrast 
to Figs. 13141 these curves are not zero-GVM contours.) 
As expected a is seen to become larger as A is reduced 
as well as when D is increased, and we found the scaling 
a oc D/A = d/ K^. a peaks when A takes values around 
the chosen fundamental wavelength Ai, and drops for 
A < Ai because the fundamental mode has maximum 
core confinement at the peak [Fig.ISf2)], while it becomes 
more poorly confined when A < Ai [Fig. [5J 1,3)]. Instead, 
the SH having A2 = Ai/2 can stay confined longer as A 
is reduced, resulting in a poor modal overlap. [A simi- 
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Table 1. Selected parameters for SHG with zero GVM. 
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